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Abstract 

Using a new definition of generalized divisors we prove that 
the lattice of such divisors for a given linear partial dif- 
ferential operator is modular and obtain analogues of the 
well-known theorems of the Loewy-Ore theory of factoriza- 
tion of linear ordinary differential operators. Possible ap- 
plications to factorized Grobner bases computations in the 
commutative and non-commutative cases are discussed, an 
application to finding criterions of Darboux integrability of 
nonlinear PDEs is given. 

1 Introduction 

Factorization is often used for simplification of solution pro- 
cedures for polynomials (factorized Grobner bases compu- 
tations) and linear ordinary differential operators (LODO). 
It is well-known that every (multivariate) polynomial fac- 
tors into product of irreducible polynomials (in the given 
coefficient field) in a unique way; for LODO an analogous 
result had been proved by E. Landau [ jL5[ and in a more 
precise form by A.Loewy M, any two different decom- 
positions of a given LODCXL into jaroducts of irreducible 
LODO L = Pi o ■ ■ ■ o Pk = Pi o ■ ■ ■ o P p have the same num- 
ber of factors (k — p) and the factors are pairwise similar 
(in some transposed order). Two (irreducible for simplic- 
ity) LODO L and M are called similar (or operators of the 
same type) if one can find operators A and B such that 
ord(A) = ord(B) < ord(L) = ord(M) and A o L = M o B 
(see below for more detail). The theory of factorization of 
LODO (Loewy-Ore theory) was developed in jl6|, [17|, |l9|. 
From the algorithmic point of view factorization of LODO 
was addressed for the first time in pi where an outline of 
an algorithm for factorization of LODO with coefficients in 
the simplest differential field of rational functions (i.e. over 
Q(a;)) was given. In the past decade many improvements of 
this algorithm and alternative algorithms were proposed (see 
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|27| ] and references therein), applications to the differential 
Galois group computation were given in [ p3[ . 

Unfortunately very little is known about factorization 
properties of linear partial differential operators (LPDO). 
The following interesting example was given by E. Landau 
(see i): if 



P = D x + xD y , 



= D X + 1, 



R = D 2 X + xD x D y + D x + (2 + x)D y , 



(1) 



then L = Q o Q o P = R o Q. On the other hand the 
operator R is absolutely irreducible, i.e. one can not factor 
it into product of first-order operators with coefficients in 
any extension of Q(x, y). 

This example shows that in order to develop a "good" 
theory of factorization of LPDO one shall try to use some 
generalization of the notion of a factor (divisor) for LPDO. 
Such tricks are very common in the commutative case, as 
the first example we may cite Kummer-Dedekind theory of 
divisors for algebraic number rings. As proposed by Dede- 
kind we may use ideals of the ring of algebraic integers of a 
given (finite) extension of Q. Since not all ideals in this ring 
are principal ideals (i.e. they are not generated as multiples 
of a single element) we obtain an extension of the notion of 
a divisor and this suffices (cf. for example JTo| , Ch.12]) to 
obtain uniqueness of decomposition of any algebraic integer 
(i.e. of the principal ideal it generates) into product of prime 
"ideal" divisors. 

For the non-commutative ring of LODO (with coeffi- 
cients in some differential field, for simplicity we will suppose 
that coefficients belong to Q(a;) i.e. they are rational func- 
tions with arbitrary algebraic number coefficients) one can 
use the Euclid division algorithm to prove that any left or 
right ideal in this LODO ring Q(a;)[,Dj is a principal ideal; 
there are no nontrivial two-sided ideals. So there is no possi- 
bility (and necessity) of "ideal" generalization of the notion 
of divisors, the only implication of non-commutativity re- 
sults in "similarity" of factors in different factorizations of 
a given LODO. 

Another well-known "unique factorization" theorem is 
the classical Jordan-Holder theorem in the theory of finite 
groups (or finitely generated modules). 

In the first half of the XX century a common approach to 
these (and many more) cases was proposed. Let us introduce 
the obvious partial order in the set of (left) ideals: Ii < I2 
if Ii D I2; we call Ii a divisor of I2 in such a case. Then 
instead of factorizations 
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L — Li o . 



(2) 
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of an element L of the ring we will consider chains \L) > 
\Li o . . . o L k ) > \L 3 o . . . o L k ) > ... > \L k ) > = ]1) of 
corresponding (left) principal ideals. Irreducibility of factors 
corresponds to maximality of this chain, i.e. impossibility to 
insert intermediate ideals between any two elements of the 
chain. The partially ordered set M (also called poset) of 
ideals in the cited above "good" cases has the following two 
fundamental properties: 

a) for any two elements A,B £ M one can find a unique 
C = sup(A,B), i.e. such C that C > A, C > B, and MX £ 
M, (X > A,X > B) X > C. Analogously there exist 
a unique D = inf(A,B), D < A, D < B, VX G M, {X < 
A,X < B) X < D. Such posets are called lattices. 
sup(A, B) and inf(A, B) correspond to the least common 
multiple and the greatest common divisor for the cases of 
number rings and LODO. In our cases a lattice will always 
have zero i.e. an element £ M such that VA" £ M, X > 
0. For simplicity (and following the established tradition) 
sup(A, B) will be hereafter denoted as A + B and inf (A, B) 

- as A ■ B; 

b) For any three A,B,C £ M the following modular 
identity holds: 



(A-C + B)-C = A-C + B-C 



(3) 



This weaker form of distributivity was discovered by De- 
dekind. The theory of modular lattices (i.e. posets with 
the above two properties, such posets are also called "Dede- 
kind structures") has beautiful (for our purpose :-) results. 
Namely these two simple properties are sufficient to prove 
the following four elegant theorems (cf. j^, [j], [LlJ]): 

Theorem 1 (Jordan-Holder-Dedekind chain condition) Any 
two finite maximal chains 

L > Li > ■ ■ ■ > L k > (4) 

L > Mi > ■ ■ ■ > M r > 
for a given L £ M have equal length: k = r. 

Theorem 2 (Kurosh & Ore) 7/ L = Li + L 2 + . . . + L p = 
M1+M2 + . . .+M r are two noncancellable sup-representations 
of L £ M then p = r and for every Li one can find Mj such 
that L = Li + ... + + Mj + L i+1 + ... + L P . 



We recall that a sup-representation 

L = Li + L 2 + . . . + L k 



(3) 



is called non-cancellable if Vi, (Li + . . . + + Li+i + 
. . . + Lk) 7^ L and each Li can not be sup-represented: Li 7^ 
Ai + Bi for Ai ^ Li, Bi ^ Li. 

We call a sup-representation (a) a direct sum if the set 
{Li} is independent that is Vi, (Li + . . . + + Li+i + 
. . . + Lk) ■ Li =0 and each Li can not be represented as a 
sum of two independent elements. Direct sums are denoted 
L = Li © L2 © . . . ffi Lk- An element A £ M is called 
indecomposable if A B @ C, B / 0, C + 0. 

Theorem 3 (O.Ore) Let an element L of a modular lattice 
have finite maximal chains fy) and L — L\ © . . . © Lj, = 
Mi © ... © M r with indecomposable Li, Mj . Then k — r and 
WLi one can find Mj such that L = L\ © . . . © Li-i © Mj © 
£;+iffi...© V 

Let us call l(L) := k + 1 the length of L £ M if L has a 
finite maximal chain mft (of length k + 1). We set 1(0) — 0. 
The length of a LODO is egual to the number of irreducible 
factors in decomposition (Q). 



Theorem 4 // all elements of a modular M have finite 
length then l(A + B) + l(A ■ B) = 1(A) + 1(B). 

These theorems give a unified approach to many well known 
facts in the theory of groups (and group representations), 
commutative and non-commutative rings; in particular they 
encompass many results of the Loewy-Ore theory of factor- 
ization of LODO. 

Let us prove here for completeness that the poset of 
(left) ideals of a (non-commutative) ring is a modular lat- 
tice. Firstly we notice that sup(j4, B) = A + B corresponds 
to the intersection of ideals A,B; inf(j4,_B) — A ■ B corre- 
sponds to the ideal composed of the sums a + b, a £ A, 
b £ B. Then if x £ A - C + B ■ C (the r.h.s of (§))) then 
x £ C, x = a + b, a £ A, a £ C, b £ B, b £ C. Obviously 
a + b£ C, a + b£ A-C+B, so A-C+B-C c (A-C+B)-C. 
Vice versa if x £ (A ■ C + B) ■ C then x £ C, x = a + b, 
a£ A, a£C,b£B=>b = x- a£C so b£B-C and 
x £ A-C + B-C which proves (A-C + B) C C A-C + B-C. 

The basic notion of similarity also exists for modular 
lattices: 

Definition 1 Two elements A, B of a modular lattice M 
are called similar if one can find C £ M such that A ■ C = 
B C = and A + C = B + C (i.e. A@C = B ®C). 

We will need also the notion of similarity of intervals or 
quotients [B/A] := {X £ M\A < X < B} for pairs A < B. 

Proposition 1 If A, B are elements of a modular lattice 
M then the intervals h = [A/(A-B)} and h = \(A + B)/B] 
are projective, I\ ~ I2, i.e. isomorphic with specific poset 
isomorphisms <f> : I\ — > I2, <t>(X) = X + B, ip(X) — A ■ X, 
iP = cf>- 1 . 

Definition 2 Two intervals [Bi/Ai], [B2/A2] are called sim- 
ilar if there exists a finite sequence of projective intervals 
[B1/A1] ~ h ~ h ~ ■ ■ • ~ h ~ [B2/A2]. 

One can prove that in Theorems the corresponding fac- 
tors (intervals) are similar (in some transposed order). Sim- 
ilarity of intervals in Theoremhl gives similarity of the re- 
spective irreducible factors in (|2]j or isomorphism of the fac- 
torgroups (factormodules) for the modular lattice of normal 
subgroups of a given (finite) group (resp. submodules). 

The case of the ring of LPDO is more complicated. It 
has no two-sided ideals and left (right) ideals are no longer 
principal ideals in the general case. Certainly the poset 
of all left (right) ideals is a modular lattice. But unfor- 
tunately we can not use the above results: for a LPDO L 
we get finite chains (^) of left ideals (the ring of LPDO 
is Noetherian, see [Q] but the intervals in any chain arc 
not (as a rule) "irreducible" i.e. one can always insert in- 
termediate ideals between some of them so the length of 
chains (W) for a given L is not bounded. For example for 



Q(x)[D x ] C Q(x,y)[D x ,D y ] we can 
> \L) + 



IDS 



arbitrary LODO L £ 
take \L) > \L) + \D™ 

= Q(x,y)[D x ,D y ]. Even the simplest D x becomes 
ducible" ! Similar infinite examples exist for decompositions 



> ...\L) + \Dy) > 



are useless, 
ideals of LPDO is 



into (direct) sup-sums. So Theorems 
We conclude that the poset of all ( 
too "large" . For the commutative case of multivariate poly- 
nomials one can limit oneself to principal ideals and get the 
desired modular lattice with finite chains. Again for LPDO 
the poset of (left) principal ideals is too "small" : it does not 
form even a lattice. For example for the two operators P, Q 
in (nl) the intersection of the left principal ideals |P) P| \Q) 
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(their "LCM") is no longer principal: one can easily check 
directly that there are no second-order common left mul- 
tiples of both P, Q but we have two linearly independent 
third-order operators divisible by P, Q: 



Ai = (xD x D y + 1 



DA 



A 



1) oP = 



(x 2 D 2 y + xD x D y - (x + 1)A ~ A) o Q, 
L 32 = (Dl + 2A + l) o P = Q o Q o P = 
{pi + xD x D y + (x + 2) A + A) o Q = R o Q, 

so there is no "least" common left multiple. Analogously we 
can directly check that these Ai , L32 have only Q, P as their 
common right divisors so Ai, A2 have no "greatest" right 
common divisor. Also as the E. Landau's example shows the 
Jordan-Holder-Dedekind chain condition fails for principal 
ideals. 

Below (section ^) we define an "intermediate" poset of 
"codimension 1" left ideals which is larger than the poset of 
principal left ideals but smaller than the lattice of all left ide- 
als. This new poset of "generalized divisors" provides all the 
necessary properties: it is a modular lattice with finite max- 
imal chains (W|) for every element and finite decompositions 
into (direct) sup-sums so the basic Theorems ^-^| are appli- 
cable; any first-order LPDO is irreducible and any LODO 
L irreducible in Q(x)[D x ] remains irreducible (as LPDO) in 
our poset. This is our main result. 

For applications the most important property of our mod- 
ular lattice of generalized divisors would be certainly the 
possibility to decompose operators into sup-sums in an overde- 
termined system of LPDO 



Lif = 0, 




A/ = 0, 


(6) 


L k f = 0. 





Suppose that A 



A 1 + 



sors then since each Ai is finitely generated (see fl 
|Ai, ■ • ■ , Li Si ), we can decompose (y) into union of 



f Ai/ = 0, 

Lisj = 0, 
A/ = 0, 

I A/ = 0. 



f W = 0> 

L PS J = 0, 
W = 0, 

I L h f = 0. 



): Ai = 
systems 



(7) 



Sums of solutions of ([?]) are obviously solutions of (|^) and 
we conjecture that they span the whole space of solutions of 
(|^) . Also we need an algorithm for such sup-decompositions 
of LPDO (see section H for the discussion) . 

Substitution of (^) with (Q) is an analogue of the well- 
known factorization technique for commutative Grobner bases 
computations. This technique considerably reduces the com- 
plexity of computations in many practical cases. The overde- 
termined systems of type (^) with one or man y u nknown 
functions are typical in many applications (Jlq, |25| ) : com- 
putation of conservation laws, symmetries ana invariant so- 
lutions of systems and single nonlinear ODEs and PDEs. 
For any system fel) one may use the standard Janet-Riquier 
technique (H ||) 

of reduction of (|^) to the so called 
passive (standard, normal) form. In the case of constant 
coefficient systems (|j) (in the commutative case) this algo- 
rithm practically coincides with the Grobner algorithm (for 



total degree+weight ordering). Unfortunately the complex- 
ity of Janet-Riquier algorithm is very high even for modest 
LPDO systems. Recently one interesting contribution to 
reduction of the complexity for computation of the genus 
(roughly speaking this is the "dimension" of the solution 
space) of (H) was given in |^]. Our approach may help in 
decomposition of the solution space of (g) into "irreducible 
submanifolds" . Further possible generalizations and appli- 
cations to the commutative case are discussed in section ^. 

Another connection of our definition of factorization of 
LPDO and integrability properties of nonlinear PDEs is dis- 
cussed in section [4f the established in |lj, |lj, |26| criterion 
of Darboux integrability [[| ("explicit" integrability) of such 
nonlinear PDEs is equivalent to generalized factorization of 
the corresponding linearized equation. This gives a new in- 
sight into possible generalizations of the notion of Darboux 
integrability of higher-order nonlinear PDEs which is now 
under investigation. 

2 Divisor ideals of LPDO 

We study general LPDO 



L : 



I i\ <m 



(8) 



\i\ = ii + . . . + i n , x = (xi, . . . ,x n ), A, = d/dxi, ord(L) := 
m. For simplicity and without loss of generality we will 
suppose that the number of the independent variables is 
n — 2, x := Xi, y := X2, and the coefficients a,ij(x,y) in (^) 
are rational functions with rational coefficients, aij{x,y) £ 
Q(x,y), so L G Q{x,y)[D x ,D y ]. 

It is straightforward to check that for every finite set of 
LPDO Li, Lk one may algorithmically find all their 
common left multiples (left cm.) up to fixed order JV: take 
Mi o Li = . . . = M k o L k with ord(Mi) = N - ord(A) and 
indefinite coefficients, then we get a linear algebraic (not 
differential!) system for the coefficients of Mi\ the number of 
equations in this linear system will be less than the number 
of the unknown coefficients for sufficiently large N, so the 
set of left (right) cm. is always nonempty. 

All these and subsequent results are certainly invariant 
w.r.t. substitution of left ideal with right ideals; application 
of the usual adjoint operation will suffice for this purpose. 
We will denote the left (right) principal ideal generated by 
LPDO L with \L) (resp. (L\). 

Definition 3 The left LPDO ideallLCMQLi) , . . . , |L fe » := 
I A) P| ■ • ■ Pi I A) is called the left least common multiple of 
LPDO Li. 

This 1LCM is always non-empty and (see Introduction) not 
principal in the general case. 

Definition 4 The same ideal |A) (")•••(") \^ k ) w ^ also 
called the left greatest common divisor of Li . 

Remark. This is serious :— ) Below the reader will see that 
this is the key to the whole trick. 

Definition 5 We call two LPDO L, R a (generalized) divi- 
sor operator couple for LPDO M if there exist LPDO X ,Y , 
Q such that 

XoM = YoR, 

X o L = Y o Q. W 
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One may informally think that M = L o R and if in fact 
M factors into the product of L, R then we may choose 
X = Q = 1, Y = Lw (|). 

A divisor operator couple is called nontrivial if ord(L) > 
0, ord(i?) > and L, R are not divisible by M, i.e. L 7^ 
M o P, R ^ K o M. In this case we will say that the 
operators M and R (M and L) have nontrivial (generalized) 
right (resp. left) common divisor. 

Remark. These definitions actually say that we can 
restore the (greatest) common divisor if we can find (least) 
common multiples; for the case of LODO if Z = 1LCM(M, R) = 
X o M = Y o R and V = rLCM(X, Y)=XoL = YoQ 
then M = LoG, R = QoG, G = rGCB(M,R). This 
explains Definition m For integral domains (commutative 
rings with 1 and no zero divisors) if two elements have LCM 
(i.e. some common multiple such that this cm. divides any 
other cm.) then they automatically have GCD (i.e. some 
common divisor such that this divisor is divided by any other 
common divisor); the converse is not true in general; also 
LCM(ac, be) = c- LCM (a, b) when one of them exist, this is 
again not true for GCD (the author thanks Dr. N.N.Osipov 
who communicated to him these facts for integral domains) . 

Lemma 1 If M) holds and for some Xi , Yi we have Xi o 
M = Yi o R, then X 1 oL = Y 1 oQ. 

Proof. Let us find some left cm. of X\, X: X = X o X\ = 
Xj o X. Then XoX 1 oM=XoY 1 o_R = Z 1 oXoM = 
X\ o Y o R so (since the ring of LPDO has no zero divisors) 
IoYj = X10Y andXioYoQ = X10X0L = ~XoX 1 oL = 
X0Y10Q hence X 1 oL = Y 1 oQ. □ 

So (0) does not depend on the choice of the 1LCM(M, R) = 
X o M = Y o R. 

Lemma 2 If ||) holds and some right cm. of M , L is cho- 
sen M o X = L o F then 



MoX 
RoX-- 



LoY, 
QoY. 



(10) 



Proof. Y o QoY = XoLoY = XoMoX = YoRoX^ 
QoY = RoX. □ 
From (g) and Lemma p] we conclude that the set of oper- 
ators L forming a generalized divisor couples with fixed R, 
M is a right ideal {Lj; from ( |10[ ) we see that for fixed M, L 
operators R form a left ideal \R}. 

Definition 6 Left ideal \R} and right ideal {L\ form a (gen- 
eralized) divisor ideal couple for an operator M (we denote 
this fact as {L\M\R}) if: 

a) any R £ \R}, L £ {L\ form a divisor operator couple for 
M i.e. (||) holds; 

b) if some LPDO L forms divisor operator couples for M 
with every R £ \R} then L £ {L\; 

c) if some LPDO R forms divisor operator couples for M 
with every L £ {L\ then R £ \R}. 



Lemma 3 Let {L\M\R} be a divisor ideal couple for M. 
Then for every Mi £ \K\ we can find a unique right ideal 
{Qi| such that {Q^M^R}. 

Proof. Since Mi£ \R} then for every L £ {L\ we have the 
unique (Lemma [l]) ft such that 3Xi, Yl, 



Take some left cm. of Mi, R for R £ \R}: Zi = X Ml 0M1 = 
Ym x R, and some left cm. of Yi, Xm 1 '■ Z2 = Yj o Xm x = 
X~ Ml o Yj . Then X Ml oX 1 o M = X Ml oY 1 o M x =Y 1 o X Ml ° 
Mi = Yi o Ya/,o R, so we get X Ml o X 1 o M = Yi o Y Ml o R. 
Using Lemma hi we conclude Xm x ° Xi o L = Yi o Ymi Q 
for any L £ {L\ and the corresponding Q. Consequently 
Fi o Y Ml °Q = X Ml oX 1 oL = X Ml o Yi 0Q1 = Yi oX Ml oQ 1 . 
Cancelling Yi we get Ym x ° Q — Xm x o Q 1 and finally 



Xmi Mi = Ym x R, 
Xm 1 o Q 1 — Yui Q, 



(12) 



(13) 



XioM = YioMi, 
Xj oL = Yi oft. 



(11) 



which shows that any ft in ([y]) (it depends on L £ {L\) 
form a divisor operator couple for Mi with any R £ \R} 
so the condition a) of Definition M holds. In order to prove 
the condition c) of it for {Qi\Ah]R} we suppose (1121 to be 
true for some LPDO R and all ft obtained from (|ll|) with 
L £ {L\. Fixing Xm x , Ymd Xi, Yi and using the same 
definition of Xm x , Y Mx : Z2 = Yi o Xmi = Xm x ° Yl we get 
FioX Ml oMi =X Ml oYioMi = YioY Ml oi? = T Ml oXioM, 
Y x oX Ml 0Q1 = Yi o Y Ml oQ = X Ml o Yi 0Q1 = X Ml 0X1 oL 

so 

X Ml o Xi o M =Yi o Y Ml o R, 
X Ml o Xi o L = Yi o Ym! o ft 

that is ([]) holds (Lemma |l|) for this R and all L £ {L\ so 
i? £ \R}- We prove the condition b) of Definition H in the 
same way: if (hj) holds for some LPDO Qi and all i? G \ R\ 
(ft Xmi, Ymi depend on R) we take R — M £ \R} so (II4) 
becomes _ _ 

Xmi o Mi = Y Ml o M, 

Xmi o Qi = Ym x o ft 

In this case we see from ([ll]) that actually we may set 
('Lemma B) X Ml = Yi, Y Ml = Xi, L := ft So we have 
QllJ), (h~2|) with i? £ and we shall prove (|) for the con- 
structed L, Q, M, R. Again Z2 = Yi o Xm x = Xa/i Yl, 
X Ml oYioAfi =I Ml oXioM =FioX Ml oMi = FioYA/jOi?, 
Xmi 0X1 oL = X Ml o Yi oft = Yi oXAf! oft = Y10 Ymi °Q, 
i.e. we have ([La) again or (Lemma hi) we have (W. □ 

Lemma ^essentially says that the right parts of {L\M\R} 
are internally characterizable as some special left ideals. The 
same is true for {L\. We will call such \R} right divisor ideals 
or r.d.i. (they are left ideals of the ring of LPDO :-) and {L 
— left divisor ideals or l.d.i. (they are right ideals). Any 
principal left ideal is a r.d.i.: \Rq) = \Rq} since we obviously 
have {l\Ro\R }, that is (|) with R = P o R , M = R , 
Y = 1, X = P, Q = P o L, (P, L being arbitrary LPDO). 

A divisor ideal couple {L\M\R} is called trivial if either 
{L\ = (M| (then \R} = = |1» or \R} = |M>, {L\ = 0. 

On the other hand as we will see in the next section there 
are divisor ideals which are not principal and not every (left) 
ideal is r.d.i. 

Let us now prove that the set of r.d.i. with the natu- 
ral ordering (\Ri} > \R 2 } iff \Ri} C |i?2» forms a lat- 
tice. Namely for two r.d.i. |-Ri}, I-R2} we take their intersec- 
tion as their 1LCM: LLCMfl-Ri}, \Ra}) := sup(|iJi}, |7? 2 }) := 
\Ri}+ \Rs} ■■= |i?i}ni^2}. Then for M £ \R!}f\\R a } 
we find two corresponding l.d.i. (Lemma |3): {L1LMI-R1}, 
{L 2 \M\R 2 }. Now let us take all LPDO L such that (§) holds 
for every R £ \Ri} f] \R 2 }. The right ideal {L\ of such oper- 
ators forms the divisor ideal couple with \R} — \Ri} f] 
for M since a) and b) in Definition hold automatically 
and {L\ D {Li|U{^ 2 l so ever y operator R such that (|9|) 
holds for all L £ {L\ forms a divisor operator couple with 



4 



all L £ {ii| and L g {L 2 \ so R £ |i?i}f||i? 2 } and the 
condition c) holds. As we will see in the next section the 
constructed {L\ is in general greater than the set of all sums 
of elements of {Li| and {£2!- Obviously this {L\ plays 
the role of lGCD^Lij, {L 2 \) = {Li| ■ {L 2 \ in the poset of 
all l.d.i. The determination of \R} = rGCD(|i?i}, \R 2 }) : = 
irif ( 1 7?i } , \R 2 }) := \Ri} • \R2} is obtained in the same way: 
we now take {L\ = rLCM({Li|, {L 2 \) := filial and 
the corresponding D |7?i}IJ|7i2} is defined using (^|). 
So our lattice of r.d.i. does not form a sublattice of the lat- 
tice of all right ideals of LPDO, it changes inf ; such subsets 
are called "meet-sublattices" . 

3 Coordinatization of divisor ideals 

Any (non-commutative) ring R satisfying the so called Ore 
condition (absence of zero divisors and existence of at least 
one common multiple for every two non-zero elements) may 
be imbedded into a skew field (non-commutative ring with 

division) built with formal quotients L _1 o M — M o L , 
L, M, T,~M £ R (Q). Let us take R = Q(x,y)[D x ] and 
the corresponding skew field Q(s, y, D x ). We can form a 
ring Q(x,y, D x )[D y ] of operators of the type L = aoDy + 
a^Dy^ 1 + . . .+a n , di — (Li) oMi £ Q(x,y, D x ) if we define 
the corresponding Dy-differentiation for the coefficients: 
d^Lii^oM^/dy = (L i )~ 1 od(M i )/dy — (L i )~ 1 oL' i o(L i )~ 1 
with d(M t )/dy = d(m {x,y)D k + ... + m k {x,y))/dy = 
(dm (x,y)/dy)D k + ... + dm k (x,y)/dy; L[ £ Q{x,y)[D x ] 
is defined via D y o Li = Li o D y + L[ =>• L' 1 o D y = D y o 
LT 1 + L7 1 oL-o Lt 1 . As explained in the basic facts of 
the Loewy-Ore theory (namely the existence of the Euclid 
division algorithm, CCDs, LCMs) hold also for operators 
with coefficients in differential skew fields. Any (left) ideal 
is again a principal ideal, they form a modular lattice. 

We have the natural projections V x : Q(x,y)[D x , D y ] — > 
Q(x,y,D y )[D x ] axidV y : Q(x,y)[D x , D y ] Q(x,y, D x )[D y ]. 
For any left ideal I C Q(x,y)[D x , D y ] we obtain two prin- 
cipal ideals V X (I) = \I X ), ~P V (I) = \I y ), their generators I x , 
I y are called coordinates of 7, I x £ Q(x, y, D y )[D x ], I y £ 
Q(x,y, D x )[D y ], we also normalize them lcof.o x (Z^) = 1, 
lcofij s (7j / ) = 1. Strictly speaking T , X {I) and V y (I) are not 
ideals in the respective rings; in order to make them ide- 
als we need to multiply all their element by various L _1 , 
L £ Q(x,y)[D x ] (resp. L £ Q(x,y)[D y }). This will be im- 
plicitly done hereafter. 

First we remark that if we have {L\M\R} then for the 
coordinates 

V x (M)=L x oR x (14) 

which gives a heuristic foundation for the definitions of the 
previous section. Also obviously ^(lLCMdHi}, I-R2})) = 
\LCM(V x (\Ri}),V x (\R 2 })) which in turn gives 
E 1 JlGCD({i 1 |,{L2|)) = lGCD(P4{ii|),P4{L 2 |)) due to 
(114). Symmetrically 

P x (rGCr>(\Ri}, |ife})) = rGCD(7>4|iii}),7M|7?2})), 
V x (iLCM{{L 1 \,{L 2 \))=rLCM(V x ({L 1 \),V x ({L 2 \)). 

The following lemma plays the key role in the subsequent 
proofs. 

Lemma 4 If a r.d.i. \R} C Q(x, y)[D x , D y ] contains two 
elements A o P, B o P such that A £ Q(a;, y) [D x ], B £ 
Q(x,y)[D y ] then P £ \R}. 

Proof. Using (Q we obtain Ao P oX = Q oY, B a P o 
X = S o Y for some 5, Q. Then P oX = A' 1 o Q oY £ 



Q(x,y,D x )[D y ], BoPoX = BoA- 1 °QoY = SoY^ 
B°A- 1 °Q = S£Q{x,y)[D x ,Dy], (15) 

B £ Q(x, y)[D x \. We will prove that in such circumstances 
Q is divisible by A: Q = AoQ,Q£ Q(x,y)[D x , D y ]. With- 
out loss of generality we may suppose lcof(B) = 1 so Jl5| ) 
reads 

(D™ + b 1 (x,y)D y n ' 1 + . . . + b m (x,y)) o 

" v ' 

B 

(16) 

(C 0J Dy + Ci-D^" 1 +... + C n )=S, 

A- 1 oQ 

with d £ Q(x,y, D y ). The leading coefficient of the l.h.s. 
(in Q(x, y, D x )[D y ]) of © is Co so since S £ Q(x, y)[D x M y ], 
Co £ Q{x,y)[D x }. Then the coefficient of D™ +n_1 in @ 
will be Ci + dCo/dy + 61 Co £ Q(x,y)[D x ] so also Ci £ 
Q(x,y)[D x ]. Using induction we get C £ Q(x,y)[D x ] so 
A^ 1 o Q — Q £ Q(x, y)[D x , D y \. This gives us the possibil- 
ity to cancel A in AoX = QoY obtaining X — QoY. Then 
the operator P £ \R} since PoX=(PoQ)oFin([lo]). □ 
Remark. Actually we used only the fact lcof_o H (B) = 1. 

Corollary 1 If a r.d.i. \R} contains two elements A, B, 
A £ Q(x,y)[D x ], B £ Q(x,y)[D y ] then \R} is trivial, \R} = 

= \1) = Q(x,y)[D x ,D y \. 

This Corollary explains why r.d.i. are "codimension 1" ide- 
als: if we will take a "codimension 2" left ideal generated by 
D x , D™, (solutions of the corresponding system D x f = 0, 
D™f = are functions of 2 variables parameterized by sev- 
eral constants and not by functions of 1 variable), it is con- 
tained only in the trivial r.d.i. |1}. The set of divisor ideals 
is a bit larger than the set of principal ideals (which are obvi- 
ously "codimension 1" ideals — solutions of the correspond- 
ing system of 1 equation are parameterized by functions of 

1 variable). 

Theorem 5 Two r.d.i. \Ri}, |7?2}, coincide iff their coor- 
dinates coincide i.e. iff P x (\Ri}) = T-^d-fi^}), P y (\Ri}) = 

-Py(\R2}). 

Proof. Let some R £ Since P x (\Ri}) = V X (\R 2 }) we 

may find R £ \R 2 }), such that C o R = R, C = (Ci)" 1 ^ £ 
Q(x, y, D y ). Multiplying with C\ we get C2 ° R — Ci o R = 
R £ \R2}- Hence for some C2 £ Q(x, y)[D y ] we have C 2 oR £ 
\R 2 }. Analogous consideration w.r.t. P y give K 2 oR £ \R 2 } 
for K 2 £ Q(x,y)[D x ], thus R £ \R 2 } (Lemma |) hence 
\Ri} C \R 2 }. Symmetrically |_R 2 } C \Ri}. □ 
Remark. This is obviously not true for arbitrary ideals: 
the ideal generated by D x , D y have the same projections as 
the trivial |1). 

Corollary 2 The lattice of r.d.i. (l.d.i.) is modular. 

Proof. Since the modular identity (0) holds for projections 
due to modularity of the lattice of the (principal) left ide- 
als of Q(x,y,D x )[D y ], Q(x,y,D y )[D x ], and GCD, LCM are 
preserved after projections, we conclude P X ((A-C+B)-C) = 
V X (A-C + B-C), P y ((A ■ C + B) ■ C) =V y {A-C + B-C), 
so (A-C + B)-C = A-C + B-C. □ 
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Corollary 3 For any r.d.i. \R} the length k + 1 of a chain 
of r.d.i. \R} > \Ri} > ... > \R k } > is limited: k + 1 < 
ord Dx (\R} x ) + ord Dy (\R} y ). 

Proof. Since V x (\Ri}), "Py(\Ri}) give chains of divisors of 
V x (\R}),V y (\R})mQ(x,y,D y )[D x ] (resp. Q(x, y, D x )[D y \), 
the adjacent elements of "P^-projected chain may differ only 
in ordo^ places (resp. in ord_o H (\R} y ) places). □ 

Analogous result is true for (direct) sup-sums of r.d.i. 

Thus Theorems lll-W are applicable to the constructed 
lattice of r.d.i. (l.d.i.j of a given LPDO. 

Proposition 2 Any first order LPDO R is irreducible (i.e. 
it has no nontrivial divisor ideal couples). 

Proof. Let R — ri(x,y)D x + r2(x,y)D y + rz(x,y) but nev- 
ertheless we have some r.d.i. \Ri}, \R} > \Ri} > 0. Neces- 
sarily ordDx(\Ri}x) = 1, ord Dy (\Ri} y ) = (up to transpo- 
sition x <-> y), so there exist Ri i = Co ° D x + G\ G -Ri}, 
Ri,2 = Ki G \Ri}, a G Q(x,y)[D y ], K x G Q(x,y)[D x ]. 
Since \R} C |.Ri}, i?i, 3 = Ri.i ~ vS^R = K 2 € \Rij, 
K2 G Q(x, y)[D y ]. Due to Corollary u\ of Lemma we have 
\Ri} = = Q(x,y)[D x ,D y ]. 

Proposition 3 Any LODO M G Q(x)[D x ] irreducible in 
this ring is irreducible as an element of Q(x, y) [D x , D y ]. 

Proof. Suppose we have a nontrivial divisor ideal couple 
{L\M\R} so for some L G {L\, R G \R} we have (Q. 
Since M G {L\, M G \R} we may suppose (subtracting 
a suitable multiple of M) that OTdD^(L) = ord(P x (L)) < 
ordo x (M) = m, ordo a .(i?) < m. Forming the quotient 
skew field Q(x,y,D x ,D y ) we can write in this skew field 
X = Y oQoL" 1 ^ X o M = Y 0Q0 L- 1 o M = Y oR=> 
Q o L -1 o M = R. If we will find lLCMCP^Q), P^(L)) = 
Z = Q°L = LoQ x G Q(a;,y,L)j / )[D I ] then ord^JL) < 
ordD a .(i) < m. For some C G Q(£, y)[D y ] we get C o 
Q = QS Q(x,y)[D x ,D v ], CoL = Z G Q(z,y) 
ord Dm (L) = ord 0j; (L) < m. Then L0Q0 L' 1 oM = LoR 
O QoLoL -1 oM = QoM = LoR. Let us take lcof _d h (L o 
7?) := lcof (7>„(Z o R)) = lcofD B (L) ■ lcofo^ (i?) = C~ o Cr G 
Q(a;,y)[£) :c ], ord(C~) < m, ord(C R ) < m, lcof Dy (Q o A/) = 
lcof u (Q) o M = C~ o Af , ord(M) = m. Finally we have 
C~ o Cjj = C~ o M in Q(x, 2/)[-Da;] which contradicts to the 

fact that M is irreducible, ord(M) = m but the orders of 
the factors in the l.h.s. are < m. □ 
Finally we are able to explain the mystery of the example 
of E. Landau ([j]): due to Proposition Q the chain \Q o Q o 
P} > \Q o P} > |P} > is maximal, so the only place 
where our generalized divisors appear is the interval between 
l-R o Q} and \Q} in \R o Q} > |Q} > 0, namely we insert 
here the I = 1LCM(|Q}, \P}). In terms of projections V X (I) 
is generated by second order (w.r.t. D x ) operator L31 = 
(x 2 D 2 y + xD x D y -(x + l)D y - D x ) o Q while T X {R o Q) 
is a third-order operator. So the operator R is reducible 
in our sense (and reducible in Q(x,y,D y )[D x ]) and it has a 
nontrivial "generalized common (right) divisor" withx 2 _D^ + 
xD x D y - (x + l)D y - D x (Definition [|. 

4 Darboux integrability of nonlinear PDEs and factoriza- 
tion of linearized equations 

In the XlX-th century vast interest in finding exact solu- 
tions to partial differential equations resulted in the devel- 
opment of methods of Lagrange, Monge, Boole and Ampere. 



G. Darboux [g] generalized the method of Monge (known as 
the method of intermediate integrals) to obtain the most 
powerful method for exact integration of partial differential 
equations known in the last century. 

Recently in a series of papers jl|, O, p6j the Darboux 
method was cast into a more precise ana efficient (although 
not completely algorithmic) form. For the case of a single 
second-order nonlinear PDE of the form 

u xy = f(x,y,u,u x ,u y ) (17) 

their idea consisted in linearization of (|l7|): using substitu- 
tion u(x,y) — > u(x,y) + ev(x,y) and cancelling terms with 
e n , n > 1, we obtain a LPDE 

V X y = AV X + BVy + CV (18) 

with coefficients depending on x, y, u, u x , u y . Equations 
of the type (|l|) were studied by Laplace, who invented 
a method of transformation (called sometimes the Laplace 
cascade method) of (|18|). First af all we remark that for the 
corresponding LPDO L = D x o D y — A ■ D x — B ■ D y — C, 

L = (D x -B)o(D y -A)+H = (D y -A)o(D m -B)+K, (19) 

where H = D X A - AB - C^K = D y B -AB-C, are called 
the Laplace invariants of (llq). So if eihter H = or K = 
our second-order LPDO L factors in the "usual" sense and 
the solutions of ( |l8| ) may be found via quadratures. If both 

H, K vanish, L is a 1LCM of two first-order LPDO. If both 
H, K do not vanish, one can apply the two Laplace trans- 
formations: L — > L\, L — > usung the substitutions 

vi = (D y - A)v, v-i = (D x - B)v. (20) 

These (invertible) transformations give two new second or- 
der LPDO Li, L_i of the same form with different coeffi- 
cients iff H 7^ (resp. K 7^ 0). In the generic case one 
obtains two infinite sequences 

L —> L\ — ► L2 — > • • • , 

L — > L-\ — > L_ 2 — ► • •• . 

If one of these sequences is finite (i.e. the corresponding 
Laplace invariant vanishes and the Laplace transform can 
not be applied once more) then the final LPDO Li is trivially 
factorable. 

One shall certainly take into consideration the original 
equation ( |l7| ) performing all the computations of the Laplace 
invariants and Laplace transforms (which allows us to ex- 
press all the mixed derivatives of u via x, y, u and the non- 
mixed U X ... X , Uy...y). 

Theorem 6 ([jl], Il4|, ^6)) A second order, scalar, hyperbolic 
partial differ entiaiequation fiF\) is Darboux integrable if and 
only if both Laplace sequences are finite. 

In EL [li], |2(| this method was also generalized for the 
case of a general second-order nonlinear PDE 

^ , y , U, U X j Uy , U XX , U X y , Uyy^J 

What can be said about factorizability of the operator L 
in our generalized sense? 

Theorem 7 L = D x o D y — a(x, y)D x — b(x, y)D y — c(x, y) 
has a nontrivial generalized right divisor ideal iff one of the 
Laplace sequences is finite. 
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Theorem 8 L — D x oD y — a(x,y)D x — b(x,y)D y — c(x,y) is 
a ILCM of two generalized right divisor ideals iff both Laplace 
sequences are finite. 

The detailed proofs will be given elsewhere. 

In fact these theorems demonstrate again that the "gen- 
eralized" factorization introduced here enjoys the necessary 
natural properties: it is invariant w.r.t. the differential sub- 
stitutions ( |2(]| ) (which destroy the "trivial" factorizations 
L = (D x - B) o (D y - A) or L = {D y - A) o (D x - B)). 

5 Conclusion 

An obvious and important generalization of our definition of 
reducibility (or sup-decomposition) of single LPDO would 
be a proper definition of decomposition of systems of LPDO 
(|^). Actually a formal generalization should be formulated 
inductively; for example if a system of 2 equations has inde- 
composable first equation Lif = then we may try to find 
(generalized) divisor operator couples for the second LPDO 
Z/2 (i.e. the second equation Lif = 0) forming for M = L2 
equations (^) modulo the left principal ideal generated by 
L\. The problem of zero-divisors (actually any LPDO is 
zero divisor now since it always has a multiple which be- 
longs to the ideal is (apparently) solved using the fact 
that for non-decomposable L\ there are no "LCM-zero divi- 
sors", i.e. if M ^ \Li), R £ \L\) then some their right cm. 
MoX = RoY \L\). In fact our proofs rely only on absence 
of "LCM-zero divisors". Certainly this approach deserves 
further thorough study in another publication. Especially 
interesting is the possibility to apply such a generalization 
to the commutative case (factorized Grobner bases). 

It would be interesting to compare our definitions of de- 
composition of ideals with the known results on decomposi- 
tion of ideals in non-commutative rings with the Ore condi- 
tion (existence of at least one common multiples for every 
two elements) 

A more challenging generalization is required for treat- 
ment of overdetermined linear partial differential systems 
with several unknown functions 

The algebraic nature of the set of operators {Q} in 
jic| ) also is of interest: we can multiply Q on the left and 
on the right with arbitrary LPDO, but addition of different 
Qi is doubly "stratified" : only Qi which belong to a fixed R 
or a fixed L may be added. As we have seen in the proof of 
Lemma H each "stratum" of {Q} is actually some \R} (resp. 

An algorithm of computation of divisor ideals for a given 
LPDO would be of big practical interest. As we have ex- 
plained in Introduction an algorithm for sup-decompositions 
is much more important for applications (Theorems EL ^|). 
One possible approach for algorithmization of sup-decom- 
positions may mimic the methods of |2J]. For th is p urpose 
we have to generalize the eigenring algorithm of []24j to the 
case of skew differential fields of coefficients Q(ie, y, D x ) with 
greater (not algebraically closed) constant subfield Q(D X ). 
Another approach to reducibility testing may rely on possi- 
ble generalization of estimates of complexity of coefficients 
of factors given in |g] for commutative coefficient fields to 
the case of the ring Q(x, y, D x )[D y ]. These difficult prob- 
lems are far beyond the scope of this short communication. 

The theorems proved in section ^ may provide a basis for 
algorithmic checking of Darboux integrability of nonlinear 
PDEs (provided a suitable factorisation algorithm for cor- 
responding linearized equations with coefficients depending 
on solutions of another PDEs will be found). Also we may 



conjecture that a generalization of the Darboux integrabil- 
ity method to PDEs of higher order with arbitrary number 
of independent variables may be given: such integrability 
should be related to representation of the corresponding lin- 
earized LPDO as a ILCM of "first-order" generalised divisor 
ideals. 
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